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This paper aims to present a general idea for description of spatially finite physical objects 
with a consistent nontrivial translational-rotational dynamical structure and evolution as a 
whole, making use of the mathematical concepts and structures connected with the Frobe- 
nius integrability/nonintegrability theorems and electromagnetic strain quantities. The idea 
is based on consideration of nonintegrable subdistributions of some appropriate completely 
integrable distribution (differential system) on a manifold and then to make use of the cor- 
responding curvatures and correspondingly directed strains as measures of interaction, i.e. 
of energy-momentum exchange among the physical subsystems mathematically represented 
by the nonintegrable subdistributions. The concept of photon-like object is introduced and 
description (including lagrangian) of such objects in these terms is given. 
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1 Introduction 

t-H ■ 

At the very dawn of the 20th century Planck (Planck 1901) proposed and a little bit later 
Einstein (Einstien 1905) appropriately used the well known and widely used through the whole 
last century simple formula E = hv,h = const > 0. This formula marked the beginning of a new 
era and became a real symbol of the physical science during the following years. According to 
the Einstein's interpretation it gives the full energy E of really existing light quanta of frequency 
v = const, and in this way a new understanding of the nature of the electromagnetic field 
was introduced: the field has structure which contradicts the description given by Maxwell 
vacuum equations. After De Broglie's (De Broglie 1923) suggestion for the particle-wave nature 
of the electron obeying the same energy-frequency relation, one could read Planck's formula in 
the following way: there are physical objects in Nature the very existence of which is strongly 
connected to some periodic (with time period T = 1/v) process of intrinsic for the object nature 
and such that the Lorentz invariant product ET is equal to h. Such a reading should suggest 
that these objects do NOT admit point-like approximation since the relativity principle for free 
point particles requires straight-line uniform motion, hence, no periodicity should be allowed. 

Although the great (from pragmatic point of view) achievements of the developed theoretical 
approach, known as quantum theory, the great challenge to build an adequate description of 
individual representatives of these objects, especially of light quanta called by Lewis photons 
(Lewis 1926) is still to be appropriately met since the efforts made in this direction, we have to 
admit, still have not brought satisfactory results. Recall that Einstein in his late years recognizes 
(Speziali 1972) that "the whole fifty years of conscious brooding have not brought me nearer to 
the answer to the question "what are light quanta", and now, half a century later, theoretical 
physics still needs progress to present a satisfactory answer to the question "what is a photon". 
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We consider the corresponding theoretically directed efforts as necessary and even urgent in view 
of the growing amount of definite experimental skills in manipulation with individual photons, 
in particular, in connection with the experimental advancement in the "quantum computer" 
project. The dominating modern theoretical view on microobjects is based on the notions and 
concepts of quantum field theory (QFT) where the structure of the photon (as well as of any 
other microobject) is accounted for mainly through the so called structural function, and highly 
expensive and delicate collision experiments are planned and carried out namely in the frame of 
these concepts and methods (see the 'PHOTON' Conferences Proceedings, some recent review 
papers: Dainton 2000; Stumpf, Borne 2001; Godbole 2003; Nisius 2001). Going not in details 
we just note a special feature of this QFT approach: if the study of a microobject leads to 
conclusion that it has structure, i.e. it is not point-like, then the corresponding constituents 
of this structure are considered as point-like, so the point-likeness stays in the theory just in a 
lower level. 

In this paper we follow another approach based on the assumption that the description of 
the available (most probably NOT arbitrary) spatial structure of photon-like objects can be 
made by continuous finite/localized functions of the three space variables. The difficulties met 
in this approach consist mainly, in our view, in finding adequate enough mathematical objects 
and solving appropriate PDE. The lack of sufficiently reliable corresponding information made 
us look into the problem from as general as possible point of view on the basis of those properties 
of photon-like objects which may be considered as most undoubtedly trustful, and in some sense, 
identifying. The analysis made suggested that such a property seems to be the fact that the 
propagation of an individual photon-like object necessarily includes a straight-line translational 
uniform component, so we shall focus on this property in order to see what useful for our 
purpose suggestions could be deduced and what appropriate structures could be constructed. 
All these suggestions and structures should be the building material for a step-by-step creation 
a self- consistent system. From physical point of view this should mean that the corresponding 
properties may combine to express a dynamical harmony in the inter-existence of appropriately 
defined subsystems of a finite and time stable larger physical system. 

The plan of this paper is the following. In Sec. 2 we introduce and comment the concept 
of photon-like object. In Sec. 3 we recall some basic facts from Frobenius integrability theory, 
then we consider its possibilities to describe interaction between/among subsystems, mathemat- 
ically represented by non-integrable subdistributions of an integrable distribution, and finally 
we introduce objects and structures in correspondence with the notion for photon-like object of 
Sec. 2. In Sec. 4 we make use of these objects to define corresponding relativistic strain tensor(s) 
and other related objects, and establish important relations with the curvature properties of the 
subdistributions considered in Sec. 3. In Sec. 5 we show explicitly how the translational-rotational 
consistency could be accounted for. In Sec. 6 we consider possible lagrangian approaches giving 
appropriate equations and spatially finite solutions with photon-like properties and behavior. 
In the concluding Sec. 6 we discuss the results obtained and present some views for further 
development. 

2 The notion of photon-like object 

We begin with the notice that any notion of a physical object must unify two kinds of properties 
of the object considered: identifying or proper, and kinematical. The identifying properties stay 
unchanged throughout the existence of the object, and the kinematical properties describe those 
changes, called admissible, which do NOT lead to destruction of the object. Correspondingly, 
physics introduces two kinds of quantities, proper and kinematical, but the more important 
quantities used in theoretical physics turn out to be the dynamical quantities which, as a rule, 
are functions of the proper and kinematical ones. In view of this we introduce the following notion 
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of Photon-like object (we shall use the abbreviation "PhLO" for "Photon-like object(s)"): 

PhLO are real massless time-stable physical objects with a consistent 
translational-rotational dynamical structure. 

We give now some explanatory comments, beginning with the term real. First we emphasize 
that this term means that we consider PhLO as really existing physical objects, not as appropriate 
and helpful but imaginary (theoretical) entities. Accordingly, PhLO necessarily carry energy- 
momentum, otherwise, they could hardly be detected. Second, PhLO can undoubtedly be 
created and destroyed, so, no point-like and infinite models are reasonable: point-like objects are 
assumed to have no structure, so they can not be destroyed since there is no available structure 
to be destroyed; creation of infinite objects (e.g. plane waves) requires infinite quantity of 
energy to be transformed from one kind to another for finite time-period, which seems also 
unreasonable. Accordingly, PhLO are spatially finite and have to be modeled like such ones, 
which is the only possibility to be consistent with their " created-destroyed" nature. It seems 
hardly believable that spatially infinite and indestructible physical objects may exist at all. 
Third, "spatially finite" implies that PhLO may carry only finite values of physical (conservative 
or non-conservative) quantities. In particular, the most universal physical quantity seems to be 
the energy-momentum, so the model must allow finite integral values of energy-momentum to 
be carried by the corresponding solutions. Fourth, "spatially finite" means also that PhLO 
propagate, i.e. they do not "move" like classical particles along trajectories, therefore, partial 
differential equations should be used to describe their evolution in time. 

The term "massless" characterizes physically the way of propagation: the integral energy 
E and integral momentum p of a PhLO should satisfy the relation E = cp, where c is the speed 
of light in vacuum, and in relativistic terms this means that their integral energy-momentum 
vector must be isotropic, i.e. it must have zero module with respect to Lorentz-Minkowski 
(pseudo)metric in M 4 . If the object considered has spatial and time-stable structure, so that 
the translational velocity of every point where the corresponding field functions are different 
from zero must be equal to c, we have in fact null direction in the space-time intrinsically 
determined by the PhLO. Such a direction is formally defined by a null vector field (, ( 2 = 0. The 
integral trajectories of this vector field are isotropic (or null) straight lines. It follows that with 
every PhLO a null direction is necessarily associated, so, canonical coordinates (x 1 , x 2 , x 3 , x 4 ) = 
(x, y,z,£ = ct) on M 4 may be chosen such that in the corresponding coordinate frame ( to 
have only two non-zero components of magnitude 1: C M = (0, 0, — e, 1), where e = ±1 accounts 
for the two directions along the coordinate z. Further such a coordinate system will be called 
zeta-adapted and will be of main usage. It may be also expectable, that the corresponding 
energy-momentum tensor T^ u of the model satisfies the relation T^ V T^ V = 0, which may be 
considered as a localization of the integral isotropy condition E 2 — c 2 p 2 = 0. 

The term "translational-rotational" means that besides translational component along £, 
the propagation necessarily demonstrates some rotational (in the general sense of this concept) 
component in such a way that both components exist simultaneously and consistently. It seems 
reasonable to expect that such kind of behavior should be consistent only with some distinguished 
spatial shapes. Moreover, if the Planck relation E = hu must be respected throughout the 
evolution, the rotational component of propagation should have time-periodical nature with 
time period T = E jv = const, and one of the two possible, left or right, orientations. It seems 
reasonable also to expect periodicity in the spatial shape of PhLO, which somehow to be related 
to the time periodicity. 

The term "dynamical structure" means that the propagation is supposed to be necessarily 
accompanied by an internal energy-momentum redistribution, which may be considered in the 
model as energy-momentum exchange between (or among) some appropriately defined subsys- 
tems. It could also mean that PhLO live in a dynamical harmony with the outside world, i.e. 
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any outside directed energy-momentum flow should be accompanied by a parallel inside directed 
energy-momentum flow. 

Finally, note that if the time periodicity and the spatial periodicity should be consistently 
related somehow, the simplest such consistency would seem like this: the spatial size along the 
translational component of propagation 4l Q is equal to cT: 4Z = cT, where l Q is some finite 
positive characteristic constant of the corresponding solution. This would mean that every 
individual PhLO determines its own length/time scale. 

We are going now to formulate shortly the basic idea inside which this study will be carried 
out. 

3 Frobenius curvature and interaction 

Any physical system with a dynamical structure is characterized with some internal energy- 
momentum redistributions, i.e. energy-momentum fluxes, during evolution. Any system of 
energy-momentum fluxes (as well as fluxes of other interesting for the case physical quantities 
subject to change during evolution, but we limit ourselves just to energy-momentum fluxes here) 
can be considered mathematically as generated by some system of vector fields. A consistent 
and interelated time-stable system of energy-momentum fluxes can be considered to correspond 
directly or indirectly to an integrable distribution A of vector fields according to the principle 
local object generates integral object. It seems reasonable to assume the following geometrization 
of the concept of physical interaction: two distributions Ai and A2 on a manifold will be said 
to interact geometrically if at least one of the corresponding two curvature forms Q1/Q2 takes 
values, or generates objects taking values, respectively in A2/A1. 

The above concept of geometrical interaction is motivated by the fact that, in general, an 
integrable distribution A may contain various nonintegrable subdistributions Ai, A2, • • • which 
subdistributions may be interpreted physically as interacting subsytems. Any physical interac- 
tion between 2 subsystems is necessarily accompanied with available energy-momentum exchange 
between them, this could be understood mathematically as nonintegrability of each of the two 
subdistributions of A and could be naturally measured directly or indirectly by the correspond- 
ing curvatures. For example, if A is an integrable 3-dimensional distribution spent by the vector 
fields (Xi, X2, X3) then we may have, in general, three non-integrable, i.e. geometrically inter- 
acting, 2-dimensional subdistributions (X±, X2), (Xi, X$), (X2, X%). Finally, some interaction 
with the outside world can be described by curvatures of distributions (and their subdistribu- 
tions) in which elements from A and vector fields outside A are involved (such processes will 
not be considered in this paper). 

To make the above statements mathematically clearer we recall the Frobenius theorem on a 
manifold M n [Godbillon 1969] (further all manifolds are assumed smooth and finite dimensional 
and all objects defined on M n are also assumed smooth). If the system of vector fields A = 
[X 1 (x),X 2 (x), . . . ,X p (x)], x G M, 1 < p < n, satisfies Xi(x) A X 2 (x) A... , A X p (x) / 0, x G M 
then A is integrable iff all Lie brackets i,j = 1,2, ...,p are representable linearly 

through the very Xi, i = 1, 2, . . . , p : [Xi, Xj] = C^X^, where C^- are functions. Clearly, an easy 
way to find out if a distribution is integrable is to check if the exterior products 

[X i ,X j ]AX 1 (x)AX 2 (x) A...,AX p (x), x G M; i,j = l,2,...,p 

are identically zero. If this is not the case (which means that at least one such Lie bracket "sticks 
out" of the distribution A) then the corresponding coefficients, which are multilinear combina- 
tions of the components of the vector fields and their derivatives, represent the corresponding 
curvatures. We note finally that if two subdistributions contain at least one common vector field 
it seems naturally to expect interaction. 
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In the dual formulation of Frobenius theorem in terms of differential 1-forms (i.e. PfafF 
forms) we look for (n — p)-Pfaff forms (a 1 , a 2 , . . . , a n ~ p ), i.e. a (n — p)-codistribution A*, such 
that (a m ,Xj) = 0, and a 1 A a 2 A • • • A a n ~ p / 0, m = 1, 2, . . . , n - p, j = 1, 2, . . . ,p. Then 
the integrability of the distribution A is equivalent to the requirements 

da m A q 1 A a 2 A • • • A a n ' p = 0, m = 1, 2, . . . , (n - p), 

where d is the exterior derivative. 

Since the idea of curvature associated with, for example, an arbitrary 2-dimensional distri- 
bution (X, Y) is to find out if the Lie bracket [X, Y] has components along vector fields outside 
the 2-plane defined by (X, Y), in our case we have to evaluate the quantities {a m , [X, Y]), where 
all linearly independent 1-forms a m annihilate (X, Y) : {a m ,X) = (a m ,Y) = 0. In view of the 
formula 

da m (X,Y) = X((a m ,Y)) - Y((a m ,X)) - (a m , [X,Y]) = -(a m , [X,Y]) 

we may introduce explicitly the curvature 2-form for the distribution A(X) = (X\, . . . ,X p ). 
In fact, if A(Y) = (Yi, . . . , Y n _ p ) define a distribution which is complimentary to A(X) and 
(a m ,Y n ) = 5™, i.e. (Yi, . . . , Y ra _ p ) and (a 1 , . . . , a n ~ p ) are dual bases, then the corresponding 
curvature 2-form Qa(x) should be defined by 

^A(X) = -da m ®Y m , since Q^X^Xj) = -da m (X u Xj)Y m = (a m , [X i ,X j ])Y m , 

where it is meant here that Oa(x) is restricted to the distribution (X±, . . . ,X p ). Hence, if we 
call the distribution (Xi, . . . , X p ) horizontal, and the complimentary distribution (Yi, . . . , Y ra _ p ) 
vertical then the curvature 2-form acquires the status of vertical bundle valued 2-form. We 
see that the curvature 2-form distinguishes those couples of vector fields inside A(X) the Lie 
brackets of which define outside A(X) directed flows, and so, do not allowing to find integral 
manifold of A(X). Clearly, the supposition here for dimensional complementarity of the two 
distributions A(X) and A(Y) is not essential for the idea of curvature. 

Hence, from physical point of view, if we make use of the quantities £l^x){Xi,Xj) to build 
the components of the energy-momentum locally transferred from the system A(X) to the 
system A(Y), then, naturally, we have to make use of the quantities ^A(y)(Y m , Y n ) to build the 
components of the energy- momentum transferred from A(Y) to A(X). It deserves to note that 
it is possible a dynamical equilibrium between the two systems A(Y) and A(X) to exist: each 
system to gain as much energy- momentum as it loses, and this to take place at every space-time 
point. On the other hand, the restriction of ^a(x) = — da m (g> Y m ,m = 1, . . . ,n — p, to the 
system A(Y), i.e. the quantities ^A(x)(Xrm ^n), and the restriction of ^a(y) = — d/3* ® Xi,i = 
1, . . . ,p, (p\ Xj) = 5), P 1 A • ■ ■ Ap p + 0, (/3 m , Y) = 0, to A(X), i.e. the quantities n A (Y)PQ, Xj), 
acquire the sense of objects causing local change of the corresponding energy-momentum, i.e. 
differences between energy-momentum gains and losses. Therefore, if W(x,Y) denotes the energy- 
momentum transferred from A(X) to A(Y), Wiy,x) denotes the energy-momentum transferred 
from A(Y) to A(X), and b~Wrx) an d SWry) denote respectively the energy-momentum changes 
of the two systems A(X) and A(Y), then according to the local energy-momentum conservation 
law we can write 

5W {X ) = W {Y ,X) + W { X,Y), 5W(y") = ~{W { X,Y) + W(y,X)) = ~SW(X)- 

For the case of dynamical equilibrium we have 6W(x) = SW(y) = 0, so we obtain 

sw {x) =o, sw (Y) = o, w {Y)X) + w {xx) = o. 

As for how to build explicitly the corresponding representatives of the energy-momentum fluxes, 
probably, universal procedure can not be offered since the adequate mathematical representative 
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of the system under consideration depends strongly on the very system. If, for example, the 
mathematical representative is a differential form G, then the most simple procedure seems to 
be to "project" the curvature components Q,M X )(Xi,Xj) and Q,^(Y)(Y m ,Y n ), as well as the 
components ^A(x)(Yi,Yj) and fi^(y)(X m , X n ) on G i.e. to consider the corresponding interior 
products. In the general case, appropriate quantities constructed out of the members of the 
introduced distributions and codistributions must be found. 

Finally we note that, as we shall see further, a PhLO may be considered to represent an 
example of a system, functioning through a dynamical equilibrium between two appropriately 
defined and interacting subsystems. 

We are going now to make use of the above general consideration to find appropriate objects 
and relations in an attempt to describe PhLO's dynamical structure and evolution in these 
terms. 



4 PhLO dynamical structure in terms of Frobenius curvature 

We consider the Minkowski space-time M = (M 4 ,ry) with signature sign{rj) = (—,—,—,+) 
related to the standard global coordinates 

(x 1 , x 2 , x 3 , x 4 ) = (x,y,z,$, = ct), 

and the natural volume form uj q = y^rfldx 1 A dx 2 A dx 3 A dx A = dx A dy A dz A d£. 

We introduce the null vector field Q 2 = 0, which in the (^-adapted coordinates (throughout 
used further) is assumed to look as follows: 

d d 

« " ^a; + w ' = ±L <» 

Let's denote the corresponding to £ completely integrable 3-dimensional Pfaff system by A*(£). 
Thus, A*(£) is generated by three linearly independent 1-forms (01,02,03) which annihilate (, 
i.e. 

«i(C) = "2(C) = "3(C) = 0; aiAa 2 AQ 3 /0. 
Instead of (a\, c*2, 03) we introduce the notation (A, A*, Q and define ( by 

C = edz + d£, (2) 

Now, since ( defines 1-dimensional completely integrable Pfaff system we have the corresponding 
completely integrable distribution (A, A*,Q. We shall restrict our further study on PhLO of 
electromagnetic nature according to the following 

Definition: We shall call a PhLO electromagnetic if it satisfies the following conditions ((, } 
is the coupling between forms and vectors): 

1. (C,A) = ((,A*) = 0, 

2. the vector fields (^4,^4*) have no components along (, 

3. (A, A*) are ^-corresponding to (A, A*) respectively . 

4. (A, A*) = 0, {A, A) = (A*, A*) . 

Remark. These relations formalize knowledge from Classical electrodynamics. In fact, our 
vector fields (^4,^4*) are meant to represent the electric E and magnetic B components of a free 
time-dependent electromagnetic field, where, as is well known [Synge,1958], the translational 
propagation of the field energy- momentum along a fixed null direction with the velocity " c" is 
possible only if the two invariants l\ = B 2 — E 2 and I2 = 2E.B are zero, because only in such 
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a case the energy-momentum tensor has unique null eigen direction. So it seems naturally to 
consider this property as intrinsic for the field and to choose it as a starting point. Moreover, in 
such a case the relation {I\) 2 + (I2) 2 = is equivalent to E 2 + B 2 = 2|E x B| and this relation 
shows that this is the only case when a nonzero field momentum can not be made equal to zero 
by means of frame change. Together with the fact that the spatial direction of translational 
energy-momentum propagation is determined by E x B, this motivates to introduce the vector 
field C in this form and to assume the properties 1-4 in the above definition. 

From the above conditions it follows that in the (-adapted coordinate system we have 

a* , , 1 9 d -j- d d 

A = uax + pay, A = —epax + euay; A = —u— p — , A*=ep— en—, 

ox ay ox Oy 

where e = ±1, and (u,p) are two smooth functions on M. 

The completely integrable 3-dimensional Pfaff system (A, A* , Q contains three 2-dimensional 
subsystems: (A,A*),(A,Q and {A*,(). We have the following 

Proposition 1. The following relations hold: 

dA A A A A* = 0; dA* A A* A A = 0; 

dA A A AC = e[u(p£ - ep z ) - p(u^ - eu z )]uj ; 
dA* A A* A( = e [u(p^ - ep z ) - p(u^ - eu z )] oj q . 
Proof. Immediately verified. 

These relations say that the 2-dimensional Pfaff system (^4,^4*) is completely integrable for any 
choice of the two functions (u,p), while the two 2-dimensional Pfaff systems (A, £) and (A*,Q 
are NOT completely integrable in general, and the same curvature factor 

R = u (ps - ep z ) - p{uj: - eu z ) 

determines their nonintegrability. 

Correspondingly, the 3-dimensional completely integrable distribution (or differential system) 
A(() contains three 2-dimensional subsystems: (A, A*), (A, () and (A*X). We have the 

Proposition 2. The following relations hold (recall that [X, Y] denotes the Lie bracket): 

[A, A*] A A A A* = 0, (3) 

d d 
[A, C] = (u ( - eu z )— + (p c - e Pz ) — , (4) 

- - d d 

= -e(p € - £ Pz)q^ + e(«£ - £u z)q^- ( 5 ) 

Proof. Immediately verified. 

From these last relations (3)-(5) it follows that the distribution (A, A*) is integrable, and it 
can be easily shown that the two distributions (A, Q and (A*, Q would be completely integrable 
only if the same curvature factor 

R = u (pi. - ep z ) - p(u/: - eu z ) (6) 

is zero (the elementary proof is omitted). 
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As it should be, the two projections 



(A, [A*,(]) = -{A*, [A, (]) = eu(ps - ep z ) - ep{u^ - eu z ) = -eK 

are nonzero and give (up to a sign) the same factor R. The same curvature factor appears, of 
course, as coefficient in the exterior products [A* , (] A A* A C and [A, (] A A A(. In fact, we obtain 

^C1a^aC = -[AC]aAac = - £ R^a|a| + rA a | a |. 

On the other hand, for the other two projections we obtain 

(A, [A,(]} = (A*, [A*,(]) = \[(u 2 +p\ - e(u 2 +p 2 ) z ]. (7) 
Clearly, the last relation (7) may be put in terms of the Lie derivative as 

l -L- c {u 2 +p 2 ) = ~L- C {A,A) = -(A,L^A) = -(A*,L;A*). 



Remark. Further in the paper we shall denote \/u 2 + p 2 = (ft>, and shall assume that ^ is a 
spatially finite function, so, u and p must also be spatially finite. 

Proposition 3. There is a function if>(u,p) such, that 

u(p£ - ep z ) - p(u£ - eu z ) _ R 



4> 2 4> 



,2' 



Proof. It is immediately verified that ij) = arctan 2 is such one. 

We note that the function ip has a natural interpretation of phase because of the easily 
verified now relations u = 4>cost^, p = (ft simp, and (ft acquires the status of amplitude. Since 
the transformation (u,p) — ► (4>,ip) is non-degenerate this allows to work with the two functions 
(4>,ift) instead of (u,p). 

From Prop. 3 we have 

K = (ft 2 L^= (ft 2 (^- £ ^ z ). (8) 

This last formula (8) shows something very important: at any (ft ^ the curvature R will 
NOT be zero only if L^ift ^ 0, which admits in principle availability of rotation. In fact, lack 
of rotation would mean that (ft and ift are running plane waves along (. The relation L^ift / 
means, however, that rotational properties are possible in general, and some of these properties 
are carried by the phase ift. It follows that in such a case the translational component of 
propagation along £ (which is supposed to be available) must be determined essentially, and 
most probably entirely, by (ft. In particular, we could expect the relation L^(ft = to hold, and if 
this happens, then the rotational component of propagation will be represented entirely by the 
phase if>, and, more specially, by the curvature factor R^O, so, since the objects we are going to 
describe have consistent translational-rotational dynamical structure, further we assume that, 
in general, L^ip ^ 0. 

We are going now to represent some relations, analogical to the energy-momentum relations 
in classical electrodynamics, determined by some 2-form F, in terms of the Frobenius curvatures 
given above. 

The two nonintegrable Pfaff systems (A,Q and (A*,^) define naturally corresponding 2- 
forms: 

G = A A C and G* = A* A (. 
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We have also the 2-vectors 



G = A A C, and G* = A* A (. 



Making use now of the Hodge *-operator, defined by r\ through the relation a A (3 = r/(*a, P)to , 
where a and (3 are p and (4 — p)-forms on M, we can verify the relation: G* = *G. The 2-forms 
G and G* define the 2-tensor, called stress-energy- momentum tensor T" according to the rule 

t ; = -\[G^G va + (GV(G*H, 
and the divergence of this tensor field can be represented in the form 

V v T;=[i{G)dG]^+[i{G*)dG*]^ 

where G and G* coincide with the metric-corresponding contravarint tensor fields, and i(G) = 
*(C) ° i(G*) = i(C) o i(^4*), i(-X") is the standard insertion operator in the exterior algebra 
of differential forms on M 4 defined by the vector field X. So, we shall need the quantities 

i(G)dG, i(G*)dG*, i(G*)dG, i(G)dG*. 

Having in view the explicit expressions for A, A*, £, A, A* and £ we obtain 



Also, we obtain 



i(G)dG = i(G*)dG* = -Lt{0 2 ).(. 



i(G*)dG = -i(G)dG* = 



(9) 



u(p£ - ep z ) - p(u^ - eu z 



dz + e 



u{Pi - £Pz) ~P(U£ - £U Z ) 



d£ = eRC- 



(10) 



In the following formulae we must keep in mind the relations d£ = 0,(^4,^4*) = ((,A*) = 
(C /!}-<). 

The two distributions (A,() and (A*,Q determine corresponding curvature forms O and O* 
according to 



n = -^dA*®A*, Q* = ~dA®A. 



Denoting Z n = tt(A,(), = Q(A*,Q, Zn* = V*(A,() and Z£„ = Sl*(A*,() we obtain 



eR 



A 



A* 



A 



eR 



(11) 



The following relations express the connection between the curvatures and the energy-momentum 
characteristics. 



i(Z n )(A A C) = 0, i(Z a )(A* A C) = eR.( = -i(G)dG* = i(G*)dG, 
i(Zn*)(A* A C) = 0, i(Z{i*)(A A £) = eR.( = -i(G)dG* = i(G*)dG, 

i(Z* n )(A A C) = 0, i(Z*n)(A* A C) = -\l^ 2 ).( = -i(G)dG = -i(G*)dG*, 
i(Z^)(A* A C) = 0, i(Z&)(A A C) = -\l^ 2 ).( = -i(G)dG = -i(G*)dG*. 
It follows from these relations that in case of dynamical equilibrium we shall have 



(12) 
(13) 

(14) 
(15) 



z(G)dG = 0, i(G*)dG* = 0, i(G*)dG + i(G)dG* = 0. 
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Resuming, we can say that Probenius integrability viewpoint suggests to make use of one 
completely integrable 3-dimensional distribution (resp. Pfaff system) consisting of one isotropic 
and two space-like vector fields (resp. 1-forms), such that the corresponding 2-dimensional spa- 
tial subdistribution (^4,^4*) (resp. Pfaff system (A, A*)) defines a completely integrable system, 
and the rest two 2-dimensional subdistributions (A,£) and (A*,() (resp. Pfaff systems (A, £) 
and (A*,()) are NON-integrable in general and give the same curvature. This curvature may 
be used to build quantities, physically interpreted as energy-momentum internal exchanges be- 
tween the corresponding two subsystems (AC) and (A*,Q (resp. (A, () and (A*,Q). Moreover, 
rotational component of propagation will be available only if the curvature R is nonzero, i.e. 
only if an internal energy-momentum exchange takes place. We see that all physically important 
characteristics and relations, describing the translational and rotational components of propa- 
gation, can be expressed in terms of the corresponding Frobenius curvature. We'll see that this 
holds also for some integral characteristics of PhLO. 

5 The electromagnetic strain viewpoint 

The concept of strain is introduced in studying elastic materials subject to external forces of 
different nature: mechanical, electromagnetic, etc. In nonrelativistic continuum physics the local 
representatives of the external forces in this context are usually called stresses. Since the force 
means energy-momentum transfer leading to corresponding mutual energy-momentum change of 
the interacting objects, then according to the energy-momentum conservation law the material 
must react somehow to the external interference in accordance with its structure and reaction 
abilities. The classical strain describes mainly the abilities of the material to bear force-action 
from outside through deformation, i.e. through changing its shape, or, configuration. The term 
elastic now means that any two allowed configurations can be deformed to each other without 
appearence of holes and breakings, in particular, if the material considered has deformed from 
configuration C\ to configuration C2 it is able to return smoothly to its configuration C\. 

The general geometrical description [Marsden 1994] starts with the assumption that an elastic 
material is a continuum BcR 3 which may smoothly deform inside the space M 3 , so, it can be 
endowed with differentiable structure, i.e. having an elastic material is formally equivalent 
to have a smooth real 3-dimensional submanifold B C M 3 . The deformations are considered 
as smooth maps (mostly embeddings) ip : B — > R 3 . The spaces B and M 3 are endowed with 
riemannian metrics G and g respectively (and corresponding riemannian co-metrics G^ 1 and 
g^ 1 ), and induced isomorphisms G and g between the corresponding tangent and cotangent 
spaces . This allows to define linear map inside every tangent space of B in the following way: a 
tangent vector V € T X M, x 6 B, is sent through the differential dip of ip to (dip) x (V) € T^^M 3 , 
then by means of the isomorphism g we determine the corresponding 1-form (i.e. we "lower 
the index"), this 1-form is sent to the dual space T*B of T^B by means of the dual linear map 
(dip)* : T*,nR 3 — > T*M, and finally, we determine the corresponding tangent vector by means of 

the isomorphism G _1 (i.e. we "raise the index" correspondingly). The so obtained linear map 

C x := [G- 1 o (dip)* o g o {dtp)] X :T X M^ T X M 

(which is denoted in [Marsden 1994] by (F T F)a;), extended to the whole B, is called Caushy- 
Green deformation tensor field. Now, the combination 

:= ^[(GoC - = off o (d^) - G] x : T t B x T t B -» R 

is called Lagrangian strain tensor field. Following corresponding linearization procedure (Ch.4 
in [Marsden 1994]) defined by an appropriate vector field X on B, representing an infinitesimal 
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displacement of B, it can be shown that the linearization of E reduces to \Lx g, where Lx is 
the Lie derivative with respect to X. 

We could look at the problem also as follows. The mathematical counterparts of the al- 
lowed (including reversible) deformations are the diffeomorphisms ip of a riemannian manifold 
(M,g), and every <p(M) represents a possible configuration of the material considered. But 
some diffemorphisms do not lead to deformation (i.e. to shape changes), so, a criterion must be 
introduced to separate those diffeomorphisms which should be considered as essential. For such 
a criterion is chosen the distance change: if the distance between any two fixed points does not 
change during the action of the external force field, then we say that there is no deformation. 
Now, every essential diffeomorphism ip must transform the metric g to some new metric (f*g, 
such that g / ip*g. The naturally arising tensor field e = (ip*g — g) ^ appears as a measure of 
the physical abilities of the material to withstand external force actions. 

Since the external force is assumed to act locally and the material considered gets the cor- 
responding to the external force field final configuration in a smooth way, i.e. passing smoothly 
through a family of allowed configurations, we need a localization of the above scheme, such 
that the isometry doffeomorphisms to be eliminated. This is done by means of introducing 
1-parameter group <pt,t € [a, b] C 1 of local diffeomorphisms, so, p a (M) and (pb(M) denote 
correspondingly the initial and final configurations. Now ipt generates a family of metrics <p>\ g, 
and a corresponding family of tensors ej. According to the local analysis every local 1-parameter 
group of diffeomorphisms is generated by a vector field on M. Let the vector field X generates 
<pt- Then the quantity 

- LvQ ■= - lim flSL — - 
2 y 2no t 

i.e. one half of the Lie derivative of g along X, is called (infinitesimal) strain tensor, or defor- 
mation tensor. 

Remark. Further in the paper we shall work with Lx g, i.e. the factor 1/2 will be omitted. 

In our further study we shall call Lx g, where g = rj is the Minkowski (pseudo)metric, just 
strain tensor. We would like to note that, as far as we know, photon-like objects have not 
been considered from such a point of view. Clearly, the term "material" is not appropriate for 
PhLO because no static situations are admissible, our objects of interest are of entirely 
dynamical nature, so the corresponding relativistic strain tensors must take care of this. 

According to the previous section important vector fields in our approach to describe elec- 
tromagnetic PhLO are (, A, A*, so, we consider the corresponding three electromagnetic strain 
tensors: Lj= m L A rj; L A * n. 

Proposition 4. The following relations hold: 



L^r] = 0, {L A n)^ = D 



2U X Uy + P X U Z Uj: 

Uy + P X 2p y P z P£ 

u z p z 





-2ep x 


-e(p y + u x ) 


-£Pz 


-ep£ 


e(p y + u x ) 


2eu y 


eu z 




-£Pz 


eu z 



















Proof. Immediately verified. 

We give now some important from our viewpoint relations. 

£>(C,C) = 0*(C,C) = o, 



n 



D(() = D{C)^ = D^Cdx" = (« c - eu z )dx + - ep 2 )dy, 
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' ' (AW 
1 



( "j - ' .^ - (p € - -7'- ^ = -[A C], 



(n 2 +p 2 ) 5 -e(u 2 +p 2 ). 



2 L 



t(p ? - ep z ) -p{u£ - eu z 



1 



-eR = — e0 2 L^. 



We also have: 



D*(Q 



(Pt ~ ep z )dx + (u(. - eu z )dy 



d d - - 

-efo - e^) — + fa -eu z )— = [A*, (], 



Clearly, and -C*(C) are linearly independent in general: 



eR 



e^ 2 L^. 



L>(C) A D*(() = e (v,£- eu z f + (p$ - ep 2 ) 2 dx Ady = ecf) 2 (ip^ - e^ z ) 2 dx Ady ^ 0. 

Recall now that every 2-form F defines a linear map F from 1-forms to 3-forms through the 
exterior product: F(a) := a AF, where a € A 1 (M). Moreover, the Hodge *-operator, composed 
now with F, gets F(a) back to *F(a) € A 1 (M). In the previous section we introduced two 
2-forms G = A A ( and G* = A* A ( and noticed that G* = *G. We readily obtain now 

£>(C) AG = D*(() AG* = D(Q A A A ( = D*{() A A* A ( = 
u(p£ — ep z ) — p(ii£ — eu z ) dx Ady Adz — u(p^ — ep z ) — p(u£ — eu z ) dx Ady A d£ = 

■ -4> 2 L^ip (e dx Ady Adz + dx Ady A d£) = -R (e dx A dy A dz + dx A dy A d£), 
D(C) AG* = -D*(() AG = D(Q AA*A( = -D*(() A A A ( = 



l r 



(u 2 +p 2 )z - e(u 2 +p 2 ) z {dx A dy A dz + e dx A dy A d£). 



Thus, recalling relations (12)- (15), we get 



D(Q A A A( = * D*(() AA* A( = -eRC = -i(G*)dG = i(G)dG 



D(Q A A* A( 



D*(()AAA( = -L^ 2 C = i{G)dG = i(G*)dG* 



(16) 



(17) 



The above relations show various dynamical aspects of the energy-momentum redistribution 
during evolution of our PhLO. In particular, equations (16-17) clearly show that it is possible 
the translational and rotational components of the energy-momentum redistribution to be rep- 
resented in form depending on the ^-directed strains D(() and D*(£). So, the local translational 
changes of the energy-momentum carried by the two vector components G and G* of our PhLO 
are given by the two 1-forms *[_D(£) A A* A £] and *[D*(Q A A A £]) and the local rotational 
ones - by the 1-forms * [D(Q A A A (] and * [D*(Q A A* A (] . In fact, the form * [D(Q AAA(~ 
determines the strain that "leaves" the 2-plane defined by (A, £) and the form * [-D*(C) AA*A( 
determines the strain that "leaves" the 2-plane defined by (A*,Q- Since the PhLO is free, 
i.e. no energy-momentum is lost or gained, this means that the two (null-field) components 
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G and G* exchange locally equal energy-momentum quantities: i{G*)dG = — i(G)dG* . More- 
over, the easily verified relation G iia G va = {G*) llcr {G*) ucr shows that the two components G and 
G* carry the same stress-energy-momentum. Now, the local energy-momentum conservation 
law V u [G^ a G ua + (G*)^ a (G*) ua ] = requires L^cp 2 = 0, and the corresponding strain-fluxes 
become zero: *[D*(() A A A (] =0, *[D(Q A A* A (] =0. On the other hand, only dynam- 
ical relation between the energy-momentum change and strain fluxes exists, so NO analog of 
the assumed in elasticity theory generalized Hooke law, (i.e. linear relation between the stress 
tensor and the strain tensor) seems to exist. This clearly goes along with the fully dynamical 
nature of PhLO, i.e. linear relations exist between the divergence terms of our stress tensor 
\[ - G^G VCJ - (G*) fMa (G*) ua ] and the C-directed strain fluxes as given by equations (16)-(17). 



6 The trans lational- rotational consistency 



We begin this section with summarizing from physical viewpoint some of the results of the 
preceding two sections in the following 

Corollary. An electromagnetic PhLO has two subsystems, mathematically represented by 
the two 2-forms G and G* , these two subsystems carry the same energy-momentum, and they 
are in a permanent dynamical equilibrium: each one gives locally to the other as much energy- 
momentum as it gains locally from it. 

This conclusion and the considerations in the preceding two sections allow to make explicit 
the mathematical representation of the PhLO translational-rotational structure. In fact, it is 
seen that the rotational component of propagation of our PhLO is dimensionally localized in a 
2-plane, which in our consideration is parametrized by coordinates (x,y), and the translational 
component of propagation is of constant nature and evolves along (. Since our PhLO is of 
electromagnetic nature, the corresponding energy-momentum tensor should be represented by 
^[ — G^G 1 " 7 — (G*) fia (G*) ua ^ . Now, the corresponding local energy-momentum conservation 
law V u [G^G vrj + {G*) pilT {G*Y' J ] = reduces to the dynamical equation L^ 2 = L^{u 2 +p 2 ) = 0, 
which seems to be naturally accepted to represent the translational component of propagation. 

In order to come to some appropriate dynamical picture of the rotational component of 
propagation we make the following consideration. Recall the two vector fields A and A*. Since 
A A A* = —e{u 2 +p 2 ) d x Ad y ^ 0, then at all space-time points occupied by our PhLO we have 
the frame Si = (A, A*, d z , d^). On the other hand the vector fields [A, (] and L4*,£] are also 
linearly independent in general: 



[A C] A[A*,(]=£ [(U£ - £U Z 



)2 + ( ^-^ )2 ]|: A | 



and we get the new frame £2 = ([ACL CL d z , d^), which is regular only if R 7^ 0. Hence, 
the rotational component of propagation transforms the frame S x to the frame S 2 . We see 
that essentially, the 2-frame (A, A*) is transformed to the 2-frame ([A, C], [A*, (]), and these two 
2-frames are tangent to the (x, y)-plane. So we get the linear map 



([A,C],[A*X]) = {A,A*) 
Solving this system with respect to (a, /3, 7, 5) we obtain 



a (3 
7 5 



a (3 


1 


-\L^ 2 sR 


7 5 




-eR -\L^ 2 



I L^ 2 
'2 <f> 2 



1 
1 



1 
-1 
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Assuming the conservation law L^<j) 2 = 0, we obtain that the rotational component of propaga- 
tion is governed by the matrix eL^ip J, where J denotes the canonical complex structure in R 2 , 
and since tp 2 L^tp = R we conclude that the rotational component of propagation is available 
if and only if the Frobenius curvature is NOT zero: R / 0. We may also say that a consis- 
tent translational-rotational dynamical structure is available if the amplitude (f> 2 = u 2 + p 2 is a 
running wave along £ and the phase ip = arctg^ is NOT a running wave along (. 

As we noted before the local conservation law L^cj) 2 = 0, being equivalent to L^cj) = 0, gives 
one dynamical linear first order equation. This equation pays due respect to the assumption 
that our spatially finite PhLO, together with its energy density, propagates translationally with 
the constant velocity c. We need one more equation in order to specify the phase function 
ip. If we pay corresponding respect also to the rotational aspect of the PhLO nature it is 
desirable this equation to introduce and guarantee the conservative and constant character of 
this aspect of PhLO nature. Since rotation is available only if L^ip ^ 0, the simplest such 
assumption respecting the constant character of the rotational component of propagation seems 
to be Lgip = const ^ 0. Now, since the usual physical dimension of the canonical coordinates 
(x,y,z,^) is [length] and the phase ip is dimensionless, we may put L^ip = const = k/I q , where 
k = ±1 and l is a positive constant with [l a ] = [length] . Note that l a is equal to the square root 
of the relation of the volumes defined by the frames £2 and £1 : l Q = ^^01(^2) /vol(T,i)\ = 
y/\aj (T, 2 )/uj (T,i)\. 

Thus, the equation L^p = and the frame rotation [A, £] = — eA* L^ip and [A*, £] = eA L^ip, 
i.e. (A, A*, d z , d%) — > ([A, (], [A*,(],d z ,d^), give the following equations for the two functions 
(u,p): 

K 

U£-EU z = - — p, p 6 

If we now introduce the complex valued function ^ 
M 2 , the above two equations are equivalent to 

which clearly confirms once again the translational-rotational consistency in the form that no 
translation is possible without rotation, and no rotation is possible without translation, where 
the rotation is represented by the complex structure J. Since the operator J rotates to angle 
a = 7r/2, the parameter l determines the corresponding translational advancement, and k = ±1 
takes care of the left/right orientation of the rotation. Clearly, a full rotation (i.e. 27r-rotation) 
will require a 4Z -translational advancement, so, the natural time-period is T = 4l /c = 1/v, 
and 4l Q is naturally interpreted as the PhLO size along the spatial direction of translational 
propagation. 

In order to find an integral characteristic of the PhLO rotational nature in action units we 
correspondingly modify, (i.e. multiply by l /c) and consider any of the two equal Frobenius 
curvature generating 4-forms: 

— dA A A A( = — dA* A 4* A ( = - eHll> = — ecf 2 Lrf> u = en-l^L u = en^- u> . 
c c cc s c l c 

Integrating this 4-form over the 4-volume M 3 x [0, 41 „] we obtain the quantity H = ekET = ± ET, 
where E is the integral energy of the PhLO and T = 4l Q /c, which clearly is the analog of the 
Planck formula E = hu, i.e. h = ET. 



K 

-£p z = —U . 

= ul + p J , where / is the identity map in 
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7 Lagrangian formulation: a complex valued scalar field 



Consider the space R 4 = R 3 x R, where R determines the time-direction. Let the field of complex 
numbers C = (R 2 , J), Jo J = — id^p. be given a real representation as a 2-dimensional real vector 
space with basis 



1 
1 



J 



1 
-1 



Every C- valued function a on R 4 can be represented in the form a = ul + pJ = <f> cos ip I + 
(ft simp. J, where u and p are two real- valued functions, cf> = \/u 2 + p 2 and ip = arctan^, and 
the components of a with respect to this basis will be numbered by latin indices taking values 
(1,2) : c*i,i = 1,2 . We denote further J (a) = —pi + uJ = a and we shall make use of the 
introduced in the preceding section constant l Q > . Now we introduce two new 1-forms: the 
first one, denoted by k s , is the restriction of j- £ to R 3 and then extended to the whole R 4 through 
zero-components in the C-adapted coordinate system: k s = j-dz\ the second one, denoted by k^, 

to 

is the restriction of 4 C to the time-direction R and then extended to the whole R 4 in the same 
way: k^ = \ Q d^. Hence, in the ("-adapted coordinate system we obtain k s = (0, 0, e/lo, 0), and 
fc£ = (0,0,0,l/Z o ). 

Making use of our vector field (, of the inner product g in C = (R 2 , J) defined by g(a, (3) = 
\tr{a o /?*), (3* is the transposed to (3, we consider now the following lagrangian (summation 
over the repeated indices: g(a,a) = oncni): 



L = 



nl g(a, L^a) + g(a, a) — nl g(a, L^a) + g(a, a) 
dai \ _ / , -„da 



Oi, 



where k = ±1. Considering a and a as independent, the Lagrange equations read 



-on 



KloC 



dx^ 



a.;. 



Note that the first equation follows from the second one under the action with J from the 
left, hence, we have just one equation of the form nl L^a = J(a), which represents the idea 
for consistent translational-rotational propagation: the translational change L^a of the field 
a is proportional to the rotational change J (a) and the coefficient l Q gives the translational 
advancement for a rotation of tt/2. 

In terms of (f> = \J u 2 + p 2 and ip = arctan^ these equations give 



L^(f>. cos(Y>) — <f>. sin(V') (^L^i> — j-^j = ^' ^C^' s i n (V0 + 4>- cos(fa (^L^tp — ^- 
These two equations are consistent only if 

L^(f> = 0, L^tfj = — . 



0. 



(18) 



The solutions are: 



<t> = (t>{x,y,£ + ez); fa = -e-z + f(x,y,£ + ez); fa = y £ + f(x, y,£ + ez), (19) 
where / is an arbitrary function. Assuming / = const we see that 



fa 



-£yZ + const = —k k^x^ + const and fa = — £ + const = k k^x^ + const 
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are the simplest possible solutions leading to non-zero curvature. We note that the spatial 
structure of the solution defined by ipi is phase dependent while the spatial structure of the 
solution defined by ^2 is NOT phase dependent. 

We note that this lagrangian leads to the obtained in the previous section linear equations 
for the components of a, which equations admit 3d-finite solutions of the kind 

a.\ = 4>cosip; ct2 = 4>smip 

with consistent translational-rotational behavior, where <j> and ip are given above, and is a 
spatially finite function. 

It is easily seen that the lagrangian becomes zero on the solutions, and since this lagrangian 
does NOT depend on any space-time metric the corresponding Hilbert energy-momentum tensor 
is zero on the solutions. This special feature of the lagrangian requires to look for another 
procedure leading to corresponding conserved quantities. A good candidate seems to be T^ v = 
(j) 2 (^( u . In fact, we obtain (in our coordinate system) 

The first term on the right is equal to zero on the solutions and the second term is zero since 
the vector field C is autoparallel, so, V^T 7 ^ = 0. 

8 Lagrangian formulation: exterior 2-form field 

We show now that the same equations for the two functions (u,p), or (4>,tp), can be obtained 
from a lagrangian defined in terms of a 2-form. Recall that the space A 2 (IR 4 ) of 2-forms on R 4 is 
6-dimensional and denote by J the complex structure in this space defined by: J\% = —J25 = 
J734 = = 3<o2 = —<J&i = 1) and all other components of J are zero in our (-adapted 

coordinate system, i.e. the only non-zero elements are the off-diagonal components, and they 
alternatively change from (+1) (upper right angle) to (—1) (lower left angle). 

We define now a representation p of the algebra C in the algebra L A 2( R 4) of linear maps in 
the 2-forms on R 4 by the relation 

dcf 

p(a e ) = p(ul + epJ) = uZ + epj, I = ^a 2 (]r 4 )i a e € C, e = ±1. (20) 

Clearly, p(a + (3) = p{a) + p(f3), p(a.(3) = p(a) o p(f3), and if G is an arbitrary 2-form then 
p(a e ).G = uG + epJ(G). Note that here and further in the text the couple (u,p) may denote 
the complex number (ul + pJ), as well as the complex-valued function a = u(x,y, z,£)I + 

Let's go back now to our (-adapted coordinate system and consider the 2-form F = dx A( = 
dx/\{edz+d^) = edxf\dz+dxf\d^. Recalling the two 1-forms A = udx+pdy and A* = —pdx+udy 
(we omit e before p and u as it was defined in Sec. 2) we obtain 

p(a £ ).F = eudx A dz + epdy Adz + udx A d£ + pdy A d£ = A A (, 

p(J(a £ )).F = {-epl + uJ).F = 

= —pdx Adz + udy Adz — epdx Ad^ + eudy A d^ = A* A ( = J(p(a £ ).F ). 

Since p is a linear map and p(0) = 0, we get one-to-one map between the C-valued functions on 
]R 4 and a special subset of 2-forms. All such 2-forms depend on the choice of the 1-form £, while 
the dependence on dx is not essential. Also, they are isotropic: 

(A A C) A (A A C) = 0, (A A C) A J (A A Q = (A A Q A (A* A 0=0, 
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i.e. they have zero invariants. 

Moreover, every such 2-form may be considered as a linear map in A 2 (M 4 ) through the above 
correspondence: p(a £ ).F a — > p(a £ ). Since together with the zero element of A 2 (M) these 2-forms 
define a linear space this property suggests to introduce inner product in this linear space 
by the rule 



< Gl(a,b),G £ (m,n) >= -tr p[a £ (m,n)] o p[a*(a,b)] 



= am + bn. 



Hence, every such 2-form acquires a norm. 

Let now F and G be two arbitrary 2-forms. In order to define the lagrangian we consider 
the Minkowski space-time M = (M 4 ,ry) as a real manifold, where the pseudoeuclidean metric 77 
has signature (—,—,—,+), and make use of the Lie derivative with respect to the vector field 
(. Also, —(k s ) 2 = (k^) 2 = (l )~ 2 - Consider the lagrangian 

L = n (k1 L^G + F, F) - 7] (k1 L^F -G,G). 

In components in the (-adapted coordinates where ( a = const, we can write 

L = \ ( Kl °~ C ^ + Fa/3 ) F ^~\ " ° a/3 ) ° <lo = C ° nSt > K = ±h (21) 

where F a/3 = 7] a ^n l3u F^ v and G alS = rj^rj^G^. Note that this lagrangian is invariant with 
respect to (F,G) —>■ (G, —F), or to (F,G) — > (—G,F). The corresponding equations read 



. ?a dG a p ^g 9F a p 



In coordinate-free form these equations look like 

k1 L^G = -F, k1 L^F = G. 

Recall now the complex structure J and assume F = —J{G). Then, treating G and J{G) 
as independent (in fact they are linearly independent on the real manifold M), and in view of 
of the constancy of J in the (-adapted coordinates, we can write 

kI q L z J{G) = -G, k1 L^G = J(G). (22) 

Since in our coordinates J and ( have constant coefficients, clearly, and J commute, so 
that the second (first) equation is obtained by acting with J from the left on the first (second) 
equation, i.e. the above mentioned invariance with respect to the transformations (F, G) — ► 
(±Cr, ^fF) is reduced to ^7-invariance. Restricting now to 2-forms of the above defined kind and 
recalling the way how J acts, J{A A Q = A* A (, i.e. the couple (^4, Q is rotated to the couple 
(A*, Q, we naturally interpret the last equations (22) as realization of the translational-rotational 
consistency: the translational change of G along ( is proportional to the rotational change of G 
determined by J, so, roughly speaking, no (-translation ( J'-rotation) is possible without 
^-rotation ((-translation), and the J'-rotation corresponds to l Q translational advancement. 

Prom (22) it follows that on the solutions the lagrangian becomes zero. So, if we try to 
define the corresponding Hilbert energy-momentum tensor the variation of the volume element 
with respect to n is not essential. Moreover, the special quadratic dependence of L on n shows 
that the variation of L with respect to n will also become zero on the solutions. Hence, this 
is another example of the non-universality of the Hilbert method to define appropriate energy- 
momentum tensor. As for the canonical energy-momentum tensor, it is not symmetric, and its 
symmetrization is, in some extent, an arbitrary act, therefore, we shall not make use of it. 
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We continue to restrict the equations (22) onto the subset of 2-forms G of the kind G = 
p(a e ).F a . As it was mentioned all these 2-forms have zero invariants: G^G^ = G ^{3 \G))^ U = 
0, or in coordinate-free way, G A G = G A J{G) = 0. Moreover, the easily verified relations 
i(()G = i(Q)J(G) = show an intrinsic connection to (: it is the only isotropic eigen vector of 
G£ = fG w and (JG)^ = ^(JG) ua . 

Substituting G = p(a £ ).F we get the already known equations 

k1 (u£ - eu z ) = -p, nl (p£ - ep z ) = u . (23) 

Clearly, (jP'Q ® Q is the right choice for energy- momentum tensor. 

Note that the 2-form F Q = dx A £ satisfies the equation L^(f> = since <^p D = 1, and does 
NOT satisfy the equation for ip, since ipF a = 0, 2ir, 4ir, so, L^(ipF ) = 0. In view of this further 
we consider only not-constant C-valued functions. 

As we already mentioned an appropriate local representative of the rotational properties of 
these solutions appears to be any of the two Frobenius 4-forms dA A A A C and dA* A A* A (, 
multiplied by the coefficient l /c, so that integrating over the 4-region (M 3 x 4l ) we get ±ET, 
which carries integral information about the rotational properties of the solution. 

We'd like to mention also that the 3-forms i(()(dA A A A Q = i(Q(dA* A A* A (), which in 
our coordinate system look like 7 A ( with 7 = —<p 2 {L^ip) dx A dy, are closed. 

The linear character of the equations obtained sets the question if the superposition prin- 
ciple holds. In general, let the parameters k, e, l Q of the two solutions be different. Let now 
Fi(k/l,£i, ll; u,p) and ^2(^2, £2, 1%; m, n) be two solutions along the same direction defined by (, 
and e of C, is of course equal to £1 for the first solution, and equal to £2 for the second solution. 
We ask now whether the linear combination c\F\ + C2-F2 with c\ = const, ci = const will be also 
a solution ^3(^3, £3, l^; c\u + cim, c\p + cin) along the same direction? In order this to happen 
the following equations must be consistent: 

K\E\l\L(U = —£ip, K\E\l\L(p = e±u, K2^ F^m = — £2?^, K2lgL^n = 82m 

n 3 e 3 llL^(ciu + c 2 m) = -e 3 (cip + c 2 n), K 3 e 3 llL^(cip + c 2 n) = £ 3 (ciu + c 2 m), 
where £3 is equal £1, or to £ 2 . The corresponding consistency condition looks as follows: 

cip + C2n c\u + C2m 



eiKlCl „ _i_ £2K2£2 „ E\K\C\ „, _|_ £2K2£2 _ " 



For example, the relations £3/^3/0 = £2^0 = £ i K i^o are sufficient for this superposition to be a 
solution. This means that if the two solutions propagate translationally for example from —00 
to +00, i.e. £1 = £2 = — 1, if the rotational orientations coincide, i.e. k± = K2, and if the spatial 
periodicity parameters are equal: l\ = l 2 = l„, the sum {c\F\ + C2-F2) gives a solution. In general, 
however, the combination c\F\ + C2F2 will not be a solution. 

On the other hand, we can introduce a multiplicative structure in the solutions of the kind 
p(a e ).F Q . In fact, if Fi(k\) = p(a 1 (ei)).F Q and -^2(^2) = p{oc 2 (£2)) -F are two such solutions we 
define their product F = F1.F2 by F = F1.F2 = p(a 1 .a 2 ).F a . Clearly, the amplitude of F is a 
product of the amplitudes of F± and F2: 4>f = 4>f 1 -4>f 2 an d the phase of F is the sum of the 
phases of F\ and F2: ipF = Y'fi + Y\f 2 - Now, F will be a solution only if 

Kp _ Kill + K2II 
£i-e 2 -£F, Tp - yrji ■ 

From the first of these relations it follows that in order F to be a solution, F\ and F2 must 
NOT move against each other, and then the product-solution shall follow the same translational 
direction. However, it is allowed F\ and F2 to have different rotational orientations, i.e. n\ / re 2) 
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then the product-solution will have rotational orientation Kp = sign(Kil 2 + fi^o) ■ Clearly, every 
subset of solutions with the same (e, k, l Q ) form a group with neutral element F Q = p(I).F Q , and 
FqT 1 = p(a^ 1 ).F , moreover, the multiplicative group of C acts on these solutions: (a £ ,F £ ) — > 
p{oi. e ).F £ . 

We give now an explicit solution, i.e. a procedure to construct a solution. We assume that 
the phase is given by ipi and e and k take values ±1 independently. The form of this solution 
u = 4> cos V'l ; p = sin ipi shows that the initial condition is determined entirely by the choice 
of <fi, and it suggests also to choose the initial condition (pt=o(x,y,ez) in the following way. Let 
for z = the initial condition <l>t=o(x , y , 0) be located on a disk D = D(x,y;a,b;r ) of small 
radius r Q , the center of the disk to have coordinates (a, b), and the value of (p t =o(x,y,0) to be 
proportional to the distance R(x, y, 0) between the origin of the coordinate system and the point 
(x, y, 0), so, R(x, y, 0) = \J x 2 + y 2 , and D is defined by D = {(x, y)| y^x — a) 2 + (y — 6) 2 < r c }. 
Also, let 6d be the smoothed out characteristic function of the disk D, i.e. #d = 1 everywhere 
on L> except a very thin hoop- like zone Bp C D close to the boundary of D where 9d rapidly 
goes from 1 to zero (in a smooth way), and 0d = outside D. Let also the dependence on z 
be given by be the corresponding characteristic function 9(z; 4Z ) of an interval (z, z + 4/ G ) of 
length 4l on the z-axis. If 7 is the proportionality coefficient we obtain 

4>{x, y, z, ct + ez) = j.R(x, y, O).0o.0(ci + ez; Al a ). 

We see that because of the available sine and cosine factors in the solution, the initial condition 
for the solution will occupy a helical cylinder of height 4Z , having internal radius of r Q and 
wrapped up around the z-axis. Also, its center will always be R(a, b, 0)-distant from the z-axis. 
Hence, the solution will propagate translationally along the coordinate z with the velocity c, and 
rotationally inside the corresponding infinitely long helical cylinder because of the z-dependence 
of the available periodical multiples. The curvature K and the torsion T of the screwline through 
the point (x, y,0) G D will be 

= R(x,y,0) = Kb 

R 2 {x,y,0) + 6 2 ' R 2 {x,y,0) + b 2 ' 

where b = 2l /ir. The rotational frequency v will be v = c/2irb, so we can introduce period 
T = 1/u and elementary action h = E.T, where E is the (obviously finite) integral energy of 
the solution defined as 3d- integral of the energy density 4> 2 . 

9 Discussion and Conclusion 

The two basic features of our approach are the assumptions for continuous spatially finite 
structure, and for available consistent translational-rotational dynamical space-time structure of 
PhLO. Hence, PhLO propagate, they do not move. The spatial structure is of composite nature 
and has two interrelated components represented by the 2-forms G and G* , and the propagation 
has two components of constant nature: translational and rotational. The translational compo- 
nent is along isotropic straight lines in (]R 4 , rj) with constant speed c. The rotational component 
of propagation is also of constant nature and follows the special rotational properties of PhLO's 
spatial structure. This intrinsically consistent dual nature of PhLO demonstrates itself according 
to the rule: no translation (rotation) is possible without rotation (translation). 

While the translational component of propagation is easily accounted through the (arbitrary 
chosen in general) isotropic autoparallel vector field £, the rotational component of propagation 
was introduced in our model making use of two things: the non-integrability properties of the 
induced by ( two 2-dimensional differential/Pfaff systems, and the complex structure J. This 
approach brought the following important consequences: 
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1. It automatically led to the required translational-rotational consistency. 

2. The rotational properties of the solutions are intrinsic for the PhLO nature, they are 
transversal to (, they are in accordance with the action of the complex structure J on the 
2-forms G and G* , and they allow the characteristics amplitude <fi and phase of a solution 
to be correctly introduced. Moreover, the rotation is of periodical nature, and a helical spatial 
structure along the spatial direction of propagation is allowed, so, the corresponding rotational 
properties differ from those in the case of rotation of a solid as a whole around a point or axis. 

3. The spatial shape and translational properties of a solution are carried by the amplitude 
(j) : L^(j) = 0, and the rotational ones are carried by the phase ip: L^ip ^ 0. 

4. The curvature, considered as a measure of the available non-integrability, or equivalently 
as a (^4* , C)-directed strain D(A*,() (or as a (A, C)-directed strain D*(A,Q) , is non-zero only 
if the phase ip is NOT a running wave along ( : L^ip ^ 0, hence, roughly speaking, curvature 
means rotation and vice versa. 

5. Quantitatively, the curvature R is obtained to be proportional to the product of the 
energy-density <p 2 and the phase change L^ip = const along ( of a solution : R = l 4> 2 . We 
recall that in General Relativity (GR) a proportionality relation between the energy-momentum 
density of non- gravitational fields and the correspondingly contracted riemannian curvature is 
postulated, while here it is obtained in the most general sense of the concept of curvature, 
namely, as a measure of Frobenius non-integrability (of correspondingly defined subdistributions 
of an integrable distribution which is meant to represent a PhLO). So, if we believe that the 
gravitational field carries energy-momentum, this result could be considered as a strong support 
for the GR ideology if the gravitational energy-momentum is manifestly included and treated in 
the same way as the energy-momentum characteristics of all other fields. 

6. A natural integral measure of the rotational properties of a solution appears to be the 
product ET, i.e. the action for one period T = 4Z G /c, which is in accordance with the Planck 
formula ET = h. 

Together with the allowed finite nature of the solutions these properties suggest the following 
understanding of the PhLO's time-stable dynamical structure: PhLO MUST always propa- 
gate in a translational-rotational manner as fast as needed in order to "survive", i.e. 
to overcome the instability (the destroying tendencies), generated by the available 
non-integrability. In other words, every free PhLO has to be able to supply immediately 
itself with those existence needs that are constantly put under the non-integrability destroying 
influence. Some initial steps to understand quantitatively this "smart" nature of PhLO in the 
terms used in the paper could be the following. 

Recalling the two 2-forms G = A A ( and G* = A* A ( we see that when L^4> 2 = then the 
corresponding subsystems keep the energy-momentum carried by each of them: i(G)dG = 
i{G*)dG* = 0. On the other hand the relation i(G)dG* = -i(G*)dG = -eR.( may be 
physically interpreted in two ways. FIRST, differentially, G transfers to G* so much energy- 
momentum as G* transfers back to G, which goes along with the previous relations stating 
that G and G* keep their energy- momentum densities. Each of these two quantities i(G)dG* 
and i(G*)dG is equal (up to a sign) to R.£, so, such mutual exchange of energy-momentum is 
possible only if the non-integrability of each of the two Pfaff 2-dimensional systems (A, Q and 
(A*, C) is present, i.e. when the curvature R is NOT zero and is measured by the same non-zero 
quantity. 

Since the curvature implies outside directed flow (with respect to the corresponding 2- 
dimensional distribution) this suggests the SECOND interpretation: the energy-momentum 
that PhLO might lose differentially in whatever way by means of G is differentially and simul- 
taneously supplied by means of G*, and vice versa. We could say that every PhLO has two 
functioning subsystems, G and G*, such, that the energy-momentum loss through the subsys- 
tem G generated by the nonintegrability of (A, (), is gained (or supplied) back by the subsystem 
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G* , and vice versa, and in doing this PhLO make use of the corresponding rotational com- 
ponent of propagation supported by appropriate spatial structure. All this is mathematically 
guaranteed by the isotropic character of G and G*, i.e. by the zero values of the two invariants 
G^ U G^ V = G*^ U G^ U = 0, and by making use of the complex structure J" as a rotation generating 
operator. 

Let's now try to express this dually consistent dynamical nature of PhLO by one object 
which satisfies one relation. We are going to consider G and G* as two vector components of 
a (M 2 , J)-valued 2-form, namely, Q = G(g)I + G*®J. Applying the exterior derivative we get 
dS7 = dG <8> I + dG* CS> J- Consider now the (R 2 , J)-valued 2- vector Q = G®I + G*®J. The aim 
we pursue will be achieved through defining the object (V is the symmetrized tensor product) 



(V, dQ) = i(G)dG ® / V I + i(G*)dG* <g> J V J + i(G)dG* + i{G*)dG 



/V J 



and put it equal to zero: (V,«)(fi,dfi) = 0. 

We note that this last relation (V, dQ) = represents the dynamical equations of the 
vacuum Extended Electrodynamics (Donev, Tashkova 1995, 2004). In particular, this equation 
contains all solutions to the Maxwell vacuum equations, and the solutions obtained in this paper 
are a special part of the full subset of nonlinear solutions to these nonlinear equations. 

From a general point of view it deserves to emphasize once again that the available curvature 
2-forms of nonintegrable subdistributions of an integrable one, represent a natural instrumentar- 
ium for describing and understanding the structure of the complex of internal exchange processes, 
which processes guarantee the time-stability of a composit physical object/system having dy- 
namical structure. The simple case of PhLO considered here suggests further applications of 
this approach to physical objects/systems of more complicated dynamical internal structure. 

In conclusion, spatially finite field models of PhLO can be built in terms of complex valued 
functions and in terms of isotropic 2-forms on Minkowski space-time, and these two approaches 
can be related. In these both cases substantial role play the complex structures J and J as 
rotation generating operators, carrying information in this way about the spin properties of 
PhLO. 

The PhLO's longitudinal sizes and rotational orientations can be determined by the constant 
parameter combination kI q , and their effective transversal sizes may be considered to depend on 
the energy-density carried by them. The time-stability is guaranteed, on one hand, by the inter- 
nal energy-momentum exchange between the two non-integrable 2-dimensional differential/Pfaff 
systems, and on the other hand, by a (possible) dynamical harmony with the outside world. 

The significance of the concepts of Frobenius curvature and (-directed strain appear to be 
of primary importance in our approach. 

The results presented in the paper clearly suggest that the Frobenius integrability theory 
appears to be important mathematical scheme to be used in theoretical physics when time- 
stable physical systems of composite structure, i.e. considered as built of relatively time-stable 
mutually interacting subsystems, are studied. The arizing Frobenius curvatures acquire natural 
interpretations of objects used to build the mathematical representatives of the corresponding 
energy-momentum exchanges, i.e. force fields. Interaction of the physical system with the 
outside world can also be incorporated in this scheme. 

We kindly acknowledge the support of the Bulgarian Science Research Fund through Contract 
$/1515. 
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